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ABSTRACT 

We perform axisymmetric (2D) multi-angle, multi-group neutrino radiation-hydrodynamic calculations of 
the postbounce phase of core-collapse supernovae using a genuinely 2D discrete-ordinate (S„) method. We 
follow the long-term postbounce evolution of the cores of one nonrotating and one rapidly-rotating 2O-M0 
stellar model for ~400 milliseconds from 160 ms to ~550 ms after bounce. We present a multi-D analysis 
of the multi-angle neutrino radiation fields and compare in detail with counterpart simulations carried out in 
the 2D multi-group flux-limited diffusion (MGFLD) approximation to neutrino transport. We find that 2D 
multi-angle transport is superior in capturing the global and local radiation-field variations associated with 
rotation-induced and SASI-induced aspherical hydrodynamic configurations. In the rotating model, multi- 
angle transport predicts much larger asymptotic neutrino flux asymmetries with pole to equator ratios of up to 
~2.5, while MGFLD tends to sphericize the radiation fields already in the optically semi-transparent postshock 
regions. Along the poles, the multi-angle calculation predicts a dramatic enhancement of the neutrino heating 
by up to a factor of 3, which alters the postbounce evolution and results in greater polar shock radii and an earlier 
onset of the initially rotationally weakened SASI. In the nonrotating model, differences between multi-angle 
and MGFLD calculations remain small at early times when the postshock region does not depart significantly 
from spherical symmetry. At later times, however, the growing SASI leads to large-scale asymmetries and 
the multi-angle calculation predicts up to 30% higher average integral neutrino energy deposition rates than 
MGFLD. 

Subject headings: Hydrodynamics, Neutrinos, Radiative Transfer, Stars: Evolution, Stars: Neutron, Stars: 
Supernovae: General 



1. INTRODUCTION 

Four decades after the first pioneering neutrino 
radiation-hydrodynamic calculations o f stellar collapse 
(IColgate & Whitd 1 19661: lArnettI 1 19661: iLeBlanc & Wilsonl 
1 19701: IWilsonI Il97l|) . the details of the core-collapse su- 
pernova explosion mechanism remain obscure. However, 
certain essentials are clear. The collapse of the evolved 
stellar core to a protoneutron star (PNS) and its evolution to 
a compact cold neutron star provides a gigantic reservoir of 
gravitational energy, ~3x 10^^ erg, a mass-energy equivalent 
of ^0.17 M0. Any core-collapse supernova mechanism 
must tap this energy and convert the fraction needed to 
match Type-II supernova observations (~10^^ erg = 1 Bethe 
[B]) into kinetic and internal energy of the exploding stellar 
envelope. 

There is general agreement that the prompt hydrodynamic 
explosion mechanism does not work and that the bounce 
shock always stalls, falling sho rt of blowing up the star (e.g., 
lBethdfT99Ql: iJankaet aDl2007|) . and must be re-energized to 
lead to a supernova. However, there is no agreement on the 
detailed mechanism that revives and endows the shock with 
sufficient energy to make a canonical ~1 -Bethe supernova. 
For decades, the "neutrino-driven" mechanism, first proposed 
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in its direct for m by [Colgate & Whitd ([T966|), and in its de - 
layed form by IWilsonI ([T985h and iBethe & Wilsonl (Il985h . 
seemed compelling. It relies on a subtle imbalance of neu- 
trino heating and cooling that leads to a net energy deposition 
behind the stalled shock, sufficient to revive it and drive the 
explosion on a timescale of hundreds of milliseconds. While 
appealing, it has been shown to fail for regular massive stars in 
spherical symmetr y (ID) when the best neutrino physics and 
transport are used (iRampp & Janka 2000; Liebendo rfer et al.l 
I2OO1I: iLiebendorfer et al.l 12005.: .Thompson et al. .2001) . Yet, 
weak explosions may be obtained i n ID for the low- 
est mass pr ogenitors, 0-Ne-Mg cores (iKitaura et al.l l2006l: 
iBurrows et'a l. 2007a). 

It is now almost certain that the canonical explosion 
mechanism must be multi-dimensional (2D/3D) in na- 
ture. The multi-D dynamics associated with convec- 
tive overturn in the postshock region (e.g., Herant et al. 



ll994l:lBurrows et al.ll 19951: 1 Janka &Mulledl 19961: iBuras et al 
2006d) and the recen tly identified standing acc r etion shock in- 
stability (SASI, e.g., Foslizzo & Tasser 2000;'Foslizzo et al.l 
[2007; S check et al. 20 ( )8; Blo ndin et al. 2003; Burrows et al.j 
I2007cl: llwakami et aP l2008l) lead to a dwell time of ac- 
creting outer core material in the postshock region that is 
larger on average than in the ID case. This results in a 
greater neutrino energy deposition efficiency behind the shock 
and, thus, creates more favorable conditions for explosion 
jBurrows & Goshv 1993|; lJankall2001l: [Thompson et"al]l2005l: 
iMarek & Janka.2007ir ^ 

The first generation of multi-dimensional supernova cal- 
culations, still employing gray flux-limited diffusion (or yet 



simpler schemes) for neutrino transport, indeed found that 
neutrino-driven convective overturn in the region between 
the stalled shock and the PNS sufficiently increased the 
neutr ino energy deposi t ion rate to lead to a delayed explo- 
sion (Herant et al. 1994; Burrows et al. 1995:'Janka & MtilleJ 
fl996: Frver & Heger 2000: Frver & Warren 2002. 2004) . The 
more sophisticated studies that followed changed this pic- 
ture. Recent long-term axisymmetric (2D) supernova calcu- 
lations with multi-group, multi-species neutrino physics and 
transport find it difficult to explode gar den-variety mas sive 
stars via the neutrino mechanism. Buras et al.l (l2006ah re- 
port explosion only f or the l ow-mass (1 1.2 M^;^) pro genitor of 
IWooslev et aP (l2002h . while lMarek & Jankal (l2007h report the 
onset of explosion in a 15-M0 model of IWoosley & Weaved 
(Il995h . given moderately fast rotation and the use of the 
Lattim er-Swesty equation of state (EOS: iLattimer & Swestyl 
Il99ll) with a nuclear compressibility modulus Kq of 1 80 MeV, 
which is significantly softer tha n the current best ex perimen- 
tal values ( fi = 240 ± 20 MeV; IShlomo et al.l[2006h . On the 
other hand. lBruenn et al.l (12006") obtain explosions for 1 l-M© 
and 15-M0 progenitors from Woosle y & Weave r (1995) only 
when they take silicon and oxygen burning into account and 
due to a synergy between nuclear burning, the SASI, and neu- 
trino heating. 

iBurrows et al.l (I2006L l2007c|) do not obtain neutrino-driven 
explosions (except in th e case of 0-Ne-Mg cores and 
accretion-induced collapse: fPessart et"ani2006bl) . but observe 
the excitation of PNS core g-modes. In their calculations, the 
PNS core oscillations reach non-linear amplitudes and damp 
via the emission of strong sound waves that propagate through 
the postshock region and efficiently deposit energy into the 
shock, eventually leading to late explosions at ~ 1 second af- 
ter bounce. This acoustic mechanism appears to be robust 
enough t o blow up even the niost massive and e xtended pro- 
genitors (iBurrows et al.ll2007d: lOtt et al.ll2006ah . but remains 
cont roversial and needs to be confirmed by other gro ups (see, 
e.g.. lYoshida et al.ll2007l: [Weinberg & Ouataertll2008h . 

In th e context of rapid pro genit or rotation. [Burrows et al.l 
(I2007bh . iDessart et al.l (l2008h . and lDessart et al.l (l2007h (^the 
latter for the accretion-induced collapse scenario) have shown 
that energetic MHD-driven explosions may be obtained 
if field-amplification by the magneto-rotational instability 
(iBalbus & Hawleylll99lh is as efficient in the core-collapse 
context as suggested (lAkiyama et al.ll2003l) . Whether rotation 
alone and without strong magnetic fields favor s or disfavors 
a neutrino-driven explosion remains to be seen (IWalder et al.l 
120051: IDessart etal.]l2006bl: lOtt et al.l |2006b|) but rapid ro- 
tation has been shown to d amp convection ([Frver & Heged 
120001) and weaken the SASI (iBurrows et al.ll2007bh . 

1.1. Core- Collapse Supernova Theory and Neutrino 
Radiation Transport 

Neutrinos, their creation, propagation, and interactions with 
supernova matter, are of paramount importance to the core- 
collapse supernova problem. They carry away ~99% of the 
final neutron star's gravitational binding energy and ~1% of 
this energy would be sufficient to blow up the star. Depend- 
ing on progenitor characteristics that set the postbounce rate 
of mass accretion onto the PNS, a successful supernova ex- 
plosion should occur within ~1-1.5 s after bounce to match 
observational and theor etical neutron star upper m ass lim- 
its around --2-2.5 M© (iLattimer & Prakashll2QQ7l and ref- 



erences therein). Consequently, the explosion mechanism 
must deliver canonical 1-B explosions on this timescale and, 
if the explosion is neutrino-driven, the neutrino heating effi- 
ciency"^ must be on the order of 10% to yield an explosion that 
achieves an energy of 1 B within ~1 s. 

The neutrinos travelling through the postshock region in 
a postbounce supernova core are not in thermal equilibrium 
with the baryonic matter. They should ideally be treated 
with full kinetic theory, describing the neutrino distributions 
and their temporal distrib ution with the Boltzmann equation 
(iMihalas & Mihalasli 19841) . Boltzmann transport is in its most 
general form a 7-dimensional problem. The 6D neutrino 
phase space (usually split up into 3D spatial coordinates, neu- 
trino energy, and 2 angular degrees of freedom) and time. In 
addition, there are up to 6 neutrino types (3 particle species, 
and their anti-particles) to deal with. Spherically- symmetric 
Boltzmann transport schemes have b een devised and imple- 
mente d in the core-collapse context (iM ezzacappa & Br uenn 
1993a;'Messer et al.iil998i; iBurrows et al. 2000: Yam ada et al 
1999; Mezzacappa & Messer 1999; Rampp & Jankai i2002i; 
iLiebendorfer et al.ii2004:iHubeny & Burrowsii2007 ). but gen- 
eral Boltzmann transport in multiple spatial dimensions is 
computationally challenging and will remain so in the inter- 
mediate term. Hence, approximations must be made in devis- 
ing computationally tractable neutrino transport schemes for 
multi-D simulations. 

A highly sophisticated approximation that arguably comes 
close to full Boltzmann transport in the cas e of quasi- 
spheric al configurations in 2D is that presented in'Buras et al.[ 
(l2006bl) , and based on earlier work by Rampp & J ankai 
(i2002h . These authors solve equations for the zeroth 
and first angular moments of spherically- symmetric radia- 
tion fields along mu ltiple radial rays (ray-by-ray approach; 
iBurrows et al."l995') and perform a variable Eddington fac- 
tor closure (Mihalas & Mihalas 1984) via a single spherically 
symmetric Boltzmann solution on an averaged ID profile 
of the 2D hydrodynamics data. Neighboring rays are cou- 
pled to provide for limited treatment of latitudinal transport. 
Their multi-group (multi-energy and multi-neutrino species) 
scheme includes inelastic neutrino-electron scattering, aber- 
ration, gravitational redshift, and frame effects to Oiy/c). 

iLivne et aP (i2004l) implemented a genuinely 2D direct so- 
lution of a reduced Boltzmann eq uation via the method 
of discrete ordinates (S^: see, e.g., iVueh & Buchler 1 977i: 
Mezzacappa & Bruennll993ai: iAdams & Larsen.2002: .Castod 



2004 and references therein) in the code VULCAN/2D, ne- 
glecting energy redistribution and fluid- velocity dependence. 

A common, more approximate way to handle neutrino 
transport that has a long pedigree in ID core-collapse studies 
is multi-group (energy/neutrino species) non-equilibrium 
flux-limited d iffusion (MGFLD; Mihalas & Mihalas 1984 . 
Arnetti il966l iBow ers & Wilsoni il982l iBruenni il985 



MvraetaLi il987i: JMvra&Burrowsi il990i: iBaron et al.i 



1989i: iCooperstein & Baroni il992i) . FLD schemes solve a 



diffusion equation for the mean radiation intensity, the zeroth 
angular moment of the specific radiation intensity. Hence, 
they drop all local angular dependence of the radiation field, 
while, in the MGFLD case, retaining the spectral neutrino 

^ We define the heating efficiency as the ratio of the energy deposition rate 
and the summed electron-neutrino and anti-electron neutrino luminosities. 
The ji and r neutrinos and their anti-particles do not contribute much to the 
heating. 



distribution. MGFLD accurately describes the radiation field 
at high optical depth where the diffusion approximation is 
exact. In the free- streaming limit, the flux must be limited to 
maintain causality and an interpolation must be performed 
between diffusion and free- streaming regimes by an ad-hoc 
prescription (using a flux limiter). 

2D F LD schemes were p i oneere d in the core-collapse con- 
text by iLeBlanc & Wilsonl (ll97Qh and modern MGFLD im- 
pl ementations ca n be fou nd in ISwesty & Myral (l2Q06h and 
in iBurrows et al.l (l2QQ7d) . It is not a priori clear whether 
MGFLD is an accurate enough prescription to yield post- 
bounce supernova dynamics in qualitative and quantita- 
tive agreement with a more accurate multi-angle treatment. 
Since net energy deposition by neutrinos is favored only in 
the semi-transparent gain layer, the quality of a MGFLD 
sche me may s ensitively depend on the flux limiter cho- 
sen (IBurrows e t al. 2000i). The fact that 2D gray FLD 
schemes have in the past led to neutrino-driven explosions 
(iHerant et al.l 19941: I Burrow s et al.ll99llFrver & Heged200Ql: 



iFrver & Warrenll20Q 2. 2004 ), while MGFLD schemes a ppear 
not to (IWalder et al.ii2005i: iBurrows et al.ll2006L l2007cl) . em- 
phasizes the importance of a spectral treatment of neutrino 
transport. 

In ID, MGFLD and Boltzmann neutrino transport 

were compared on s tatic hydr odynamic po st bounc e 

backgrounds bv iJankal (Il992h. lYamada et al.l (Il999l) . 
iMesser et a l.' (1199 81). and IBurrows et al.1 (l200Ql) . Also in 
ID, iMezz acappa & BruennI (Il993b|) compared Boltzmann 
transport and MGFLD ev olutions in the collapse phase, while 
[Liebendorfer e t al.l (l2004l) performed the only comparison to 
date of ID long-term Boltzmann and MGFLD supernova evo- 
lutions. The static studies all agree that Boltzmann transport 
yields larger instantaneous neutrino heating rates in the gain 
region, mostly because of a more slowly decreasing inverse 
flux factor (c over the ratio of flux to neutrino energy density), 
a quantity that can be related to the rat e of en ergy absorption. 
On the other hand, iLiebendorfer et al.l (l2004l) find no signifi- 
cant dynamical differences between MGFLD and Boltzmann 
transport evolutions in their long-term comparison study 
with t he 13-M0 progenitor model of Nomoto & Hashimota 
([19881) . 

In this paper, we present 2D multi-angle, multi-group neu- 
trino transport supernova calculati ons using the Newtonian 
axisy n imetric VULCAN/2D code (iLivnd 119931: iLivne et al.l 
I2Q04I: IBurrows et aP l2007cl) . Comparing multi-D Boltz- 
mann and MGFLD treatments, we perform postbounce sim- 
ulations with VULCAN/2D and compare 2D steady- state 
snapshots, as well as fully-coupled dynamical 2D radiation- 
hydrodynamics evolutions, for non- and rapidly-rotating 
20-Mc7) models who se precollapse profiles are taken from 
IWoosley et al.l (l2002h . We analyze our angle-dependent neu- 
trino radiation fields and provide for the first time local 2D 
map projections of the specific intensity Ij^. 

In ^ we describe our hydrodynamic and radiation- 
transport schemes and the microphysics that we use in this 
postbounce core-collapse supernova study. In ^ we intro- 
duce the presupernova models and the postbounce configu- 
rations, the setup, and the methodology of our Boltzmann- 
transport-MGFLD comparisons. In ^4) we present results of 
snapshot Boltzmann transport calculations and compare them 
with their MGFLD counterparts. In ^ we then discuss time- 
dependent calculations, the dynamical differences between 



Boltzmann transport and MGFLD runs, and the consequences 
for postbounce supernova model evolution. We wrap up in ^ 
with a summary and critical discussion of the work presented 
in this paper. 

2. METHODS 

2.1. Hydrodynamics 

We employ the arbitrary Lagrangian-Eulerian (ALE, 
with second-order total- variation-diminishing [TVD] remap) 
radiation-hydrodynamics code VUL CAN/2p. The hydrody- 
namics module was first described by lLivnd (1 19931) ^. The 2D 
time-explicit hydrodynamics scheme is second-order accurate 
(in smooth parts of the flow), unsplit, and implements a finite- 
difference representation of the Newtonian Euler equations 
with artificial viscosity on arbitrarily structured grids and in 
cylindrical coordinates. The computational grid employed 
here is set up to resemble a spherical-polar grid at radii greater 
than 20 km an d gradu ally transitions to a Cartesian structure at 
smaller radii (lOtt et a l. 2004). This (a) avoids hydrodynamic 
timestep restrictions due to focussing of angular grid lines and 
(b) liberates the PNS core, thus allowing mass motion along 
the axis of symmetry. 

Self-gravity is implemented via direct grid-based solu- 
tion of the Newtonia n Poisson equation, as described in 
IBurrows et al.l (l2007d) . and \ ye employ the finite -temperature 
nuclear equation of state of IShen et al.l (Il998allbl) . The cal- 
culations are run with 230 logarithmically- spaced radial and 
120 angular zones (including the inner, quasi-Cartesian re- 
gion). The grid encompasses a radial extent of 4000 km and 
the full 180° of the axisymmetric domain. 

2.2. Neutrino Transport and Microphysics 

VULCAN/2D contains two multi-group, multi- species neu- 
trino radiation-transport options. As we discuss below (^, 
both modules are used in this study. The module imple- 
menting 2D transport in the MGFLD approximation, evolv- 
in g the zeroth moment o f the radiation field, is discussed 
in IBurrows et al.l (l2007d) . The angle-dependent transport 
module that evolves the specific neutrino radiation intensity, 
/(r, (^, £j,, species, t), via the method of d iscrete ordi nates (S^), 
was first discussed bylLivne et al.l (12004)) (se e also iMorel et al.l 
[T996I: lAdams & Larsenll2002l:ICastoJl2004l) . 

For convenience and future reference, we define the zeroth, 
first, and second moments of the radiation field. 






Hjy = — (p dVt nly , 

47r/4^ 

Kjj = — (i) dft nnly . 

47r/4^ 



(1) 
(2) 
(3) 



Note the vector and tensor natures of ^^^ and K^^, respectively, 
^is the radiation field unit vector whose coordinate-dependent 
components are given in Hubeny & Burrows (2007) for vari- 
ous common coordinate systems. Here we employ cylindri- 
cal coordinates (see Fig. [T]). The radiation-pressure tensor 

^ For details and an exte nsion to magneto-hydrodynamics not employed 
here, see lLivne et al.1 J2007D . 



Kjy obeys the trace condition Jjy = Tr[Kjy] (iMihalas & Mihalai 
Il984l) . The spectral neutrino flux is defined as /v = ^"kHv 

As explained in iLivne et aP (l2004l) . the time-implicit S„ 
solver in VULCAN/2D updates the specific intensity in 
the laborator y frame via the Boltzmann transport equation 
(ICastodll972h without fluid- velocity dependence, 



- — +n-VI+aI = S . 
c dt 



(4) 



where we have dropped the neutrino group index v. a = a^ + 
a\ where cr^Cr,^^^, species) is the inverse absorption mean- 
free path and cr^(r, e^, , species) is the inverse scattering mean- 
free path (both equivalent to the corresponding cross sec- 
tion multiplied by the number density). We assume scatter- 
ing to be isotropic and employ the transport cross section 
cr^ = (1 - {cos i9))a^j^ instead of the total scattering cross sec- 
tion a^j^. This approach has been shown to work well in 
spherically- symmetric core-collapse supern ova calculations 
(iBurrows et al. I I200QI: [Thompson et al.ll2003h . The right-hand 
side source term S equals 5'em(i',^zy, species) + cr^/, where 5'em 
is the emissivity. The transport grid is identical to the hydro- 
dynamics grid. The specific intensity and its moments are 
defined at cell centers, facilitating spatially-consistent cou- 
pl ing with the scalar hydrodynamics variables, as discussed 
in lLivne et al.l (12004)) . Radiation stress at cell corners is com- 
puted via linear interpolation employing cell-centered values 
of the radiation flux. 

As a consequence of the neglect of Oiy/c) terms in our 
transport formulation, neutrino advection, Doppler shifts and 
aberration effects are not considered. This greatly limits 
the computational complexity of the problem, but its impact 
on the transport solution depends on t he particular choice 
of refe rence frame and was examined in iHubeny & Burrowsl 
(|2007|) . It is clear that around core bounce and neutrino 
breakout, during the non-linear phase of the SASI hundreds 
of milliseconds after bounce, and in the case of rapid rota- 
tion, including 0(y/c) terms is advisable. We leave them 
out here in order to make long-term multi-angle radiation- 
hydrodynamics simulations feasible and allow direct compar- 
ison with the MGFLD variant of VULCAN/2D. Full 0(y/c) 
Boltzmann transport with energy redistribution will be ad- 
dressed using the code BETHE currently under develop- 
ment by a subset of our group (iHubeny & Burrowsl 120071: 
iMurphv & Burrowsll2008h . 

We discretize the angular radiation distribution evenly in 
cos^ from -1 to 1 and in (p evenly from to tt (treating only 
one hemisphere because of axial symmetry). We make the 
number of (/:)-bins a function of cos "d to tile the hemisphere 
more or less uniformly in solid angle. In our time-dependent 
runs we employ 8 cos^ bins, resulting in a total of 40 angular 
zones. Steady-state radiation fields are computed with 8 cos^ 
bins, 12 cos^ bins (92 total angular zones) and 16 cos^ bins 
(162 total angular zones) at each spatial grid point. 

The standard set of neutrino-matter interactions listed in 
[Thompson et al.l (l2003h is included and all computations are 
performed with 16 discrete neutrino energy bins, approxi- 
mately logarithmically spaced from 2.5 MeV to 220 MeV. 
Electron neutrinos (Ve) and electron- antineutrinos (Ue) are 
treated independently while we lump together the heavy- 
lepton /i, /i, r, and f neutrinos into one group ('V^"). The 
code is very efficiently parallelized via MPI in energy groups 
and species. As an additional simplification, we do not in- 




FiG. 1 . — Coordinates used in the axisymmetric S^ transport scheme imple- 
mented in VULCAN/2D. The radiation direction vector nis defined in terms 
of i9 and (yi?. 1? is the angle with respect to the coordinate-grid z-axis at all spa- 
tial positions (z,cn). At each (z,03), the local momentum- space unit sphere is 
covered by n zones in i? and at each i^ location by a number m(i&) of c/^-zones, 
so that each zone in (i?, vp) covers roughly the same solid angle. 
elude energy redistribution by inelastic neutrino-electron scat- 
tering. Such energy redistribution and scattering are of mod- 
est (^ 10%) relevance for the trapped electron fraction (Ye) 
and entropy of the inner co re at core bounce, but otherwise 
arguably quite subdominant ([Thompson et al.|[2003|) . 

2.3. A Hybrid Approach - Combining S„ and MGFLD 
Neutrino Transport 

The time-implicit S„ scheme in VULCAN/2D is iterative 
and suffers convergence problems in regions where the trans- 
port problem is scattering-dominat ed and the optical depth 
is high (r>5). As a consequence, iLivne et aP (l2004|) lim- 
ited the timestep at postbounce times to ~0. 1-0.3 fis to en- 
sure accuracy and stability. In the present study, we take a 
different approach and introduce a hybrid S„-MGFLD trans- 
port scheme that treats the quasi-isotropic transport problem 
in the optically-thick PNS interior in the diffusion approxima- 
tion and transitions to full multi-angle S„ transport in a region 
of moderate optical depth (r > 2), but that is still significantly 
interior to the neutrinospheres (r ~ 2/3) where the neutrinos 
decouple from matter and begin to stream. 

We chose a radius of 20 km in our calculations for the tran- 
sition from MGFLD to S„. This is a sensible choice, (1) be- 
cause the neutrinosphere radii of all groups (energies/species) 
remain larger than 20 km throughout the postbounce period 
our simulations cover and, (2) because 20 km also marks the 
radius at which the transition from the inner irregular quasi- 
Cartesian grid to the outer regular grid is complete. This 
boundary is smooth and the S„-MGFLD transition does not 
suffer from Cartesian cornerstone effects. 

The transition is implemented by setting up for each energy 
group and species an approximate specific intensity /j, at the 
centers of the zones below the S„-MGFLD interface using 
the information available from MGFLD. This approximate /^y 
is obtained via its angular expansion to first order in n (the 
Eddington approximation): 

U=Io + 3(nH). (5) 

Here, /q = /mgfld and H = #MGFLD/47r, where Fmgfld is the 
flux, and H is the first moment of /j,. In MGFLD, Fmgfld is 
computed via 

^MGFLD=-FL[Z)]V/, (6) 
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Fig. 2. — Entropy colormaps of the nonrotating model s20.nr (left) and the rotating model s20.7r (right) at 160 ms into their postbounce evolution computed 
with MGFLD. Velocity vectors are superposed with vector lengths saturated at l.Ox 10^ cm s~^ . Model s20.nr has a practically spherical PNS and shows features 
of violent overturn in the convectively unstable postshock region. The shock radius in this model is ~ 175 km at this point and the onset of the SASI is apparent 
from the slightly deformed shock. Model s20.7r, on the other hand, has a strongly rotationally-flattened PNS and convective overturn is confined to polar regions. 
These regions exhibit the globally highest entropies and greatest entropy gradients, since the polar velocity divergence at the shock is the highest. The shock 
radius at this time in model s20.7r is ~230 km and no SASI features are visible. 



where 



D = 



1 

3^ 



(7) 



(8) 



with Bruenn's flux Hmiter^ (lBrueniilll985h . 
FL[Z)]= ^ . 

The first angular moment of eq. (|5]), fs„ = j nIdVt, is then 

equar to /^mgfld and the S„-MGFLD matching is consistent 
and provides a representation of the specific intensity / that 
is accurate to first order in n. Given the essentiaUy isotropic 
neutrino radiation field deep inside the PNS, this approxima- 
tion yields excellent results. We note that the scheme makes 
the implicit assumption that the radial gradient of the mean 
intensity at the transition radius is always negative or zero. 
This condition is generally fulfilled in PNSs. 

3. INITIAL MODELS AND SETUP 

We employ the spherically- symmetric solar-metallicity 
2O-M0 (at ZA MS) model s2Q fro m the stellar evolu- 
tionary study of IWoosley et alJ (l2002h . who evolved it to 
the onset of core collapse. At that moment, its iron core 
mass^ is ~1.46 M© and its central density has reached 



^ W e use Bruenn's flux limiter in VULCAN/2D, because iBurrows et al.l 
(1200(1) found it to perform best in their comparison of flux limiters with angle- 
dependent transport. 

^ Given the limited number of angular zones of / and the fact that we are 
not using Gaussian-quadrature-type angular zoning, the integrals of / are only 
accurate to ~5% when 8 'j?-zones are used and accurate to ~1% when 12 and 
more ^-zones are employed. To ensure conservation of energy in the S^- 
MGFLD matching, we employ purely geometrical and temporally constant 
correction factors to enforce F^^ = /mgfld at the interface. 

^ Determined by the discontinuity in the electron fraction. Ye, at the outer 
edge of the iron core where 7^^0.5. 



~8.4x 10^ g cm~^. A graph of the progenitor's precollapse 
density stratificatio n as a function of enc losed mass can be 
found i n Fig. 1 of [B urrows et al. I (l2007d) . Note that in the 
study of lWoosley et al. (2002) , iron core mass and extent vary 
non-monotonically in the 10-20 M© ZAMS mass range and 
that their solar-metallicity 2O-M0 model has, in fact, a more 
compact central configuration than the corresponding IS-M© 
model. Stellar evolution theory of massive stars has yet to 
converge and studies by different groups do not presently 
yield the same presupernova structures. 

We set up two initial models in VULCAN/2D: s20.nr and 
s20.7r. Both models are mapped from ID onto our 2D hydro- 
dynamic grid under the assumption of spherical symmetry. 
Model s20.nr is kept nonrotating, while we impose an initial 
angular velocity profile in model s20.7r according to the rota- 
tion law, 

where 03 is the distance from the rotation axis and A is a pa- 
rameter governing precollapse differential rotation. This ro- 
tation law enforces constant angular velocity on cylindrical 
shells and, for se nsible choices of A, reproduces qualitatively 
(lOtt et al.ll2006b|) predictions from presupernova models tha t 
include rotation in a ID fashion (iHeger et al.l l2000L l2005h . 
Since the computational complexity of this study inhibits us 
from performing a sweep of the Qq-A parameter space, we 
chose A = 1000 km and Qq = n rad s~^. Hence, the initial 
central period is 2 s - this is an identical r otational setup to 
the fiduc i al mode l in lBurrows et al.l (l2007bh . As discussed in 
lOtt et al.l (l2006bh . 2 s is rather short, leads to a rapidly rotat- 
ing postbounce configuration with a millisecond-period PNS, 
and, unless significant postbounce spin-down (e.g. via MHD 
torques) occurs, is inconsistent with average pulsar birth spin 
estimates. We chose such rapid rotation simply because we 



TABLE 1 
Model Summary 



Model Name 


Progenitor 


(rads-i) 


A 
(km) 


tb 
(ms) 


^snap 

(ms) 


tf-tt 
(ms) 


s20.nr 
s20.7r 


s20.0 
s20.0 


0.0 

TT 


1000 


179.2 
193.7 


160.0 
160.0 


500.0 
550.0 



Note. — Summ ary of model parameters. The progenitors are taken from 
IWooslev et alj i2002|) . Qq is the initial central angular velocity, A is the differential 
rotation parameter of the rotation law (eq. [9). th is the time of core bounce, ?snap is 
the time after tb at which the postbounce snapshots are taken, and tf - tb is the point at 
which we stop our simulations. 

wish to study a postbounce supernova core with significant 
rotationally-induced asymmetry. Key model parameters and 
characteristics are summarized in Table [B 

We collapse both models with the MGFLD variant of VUL- 
CAN/2D and evolve them to ~160 ms after core bounce. 
Then, we transition to S„ Boltzmann transport and solve for 
the stationary neutrino radiation field based on the artificially 
frozen hydrodynamics data at this postbounce time. Once we 
have obtained a converged angle-dependent radiation field, 
we activate neutrino-matter coupling and hydrodynamics and 
evolve in time the coupled radiation-hydrodynamics equa- 
tions. For direct comparison, we also continue the MGFLD 
simulations to later times. All steady-state snapshots are com- 
puted in three momentum- space angular resolutions, Si6, S12, 
and Ss, while the long-term evolution calculations could only 
be performed with Sg, due to computational constraints. 

In Fig. [21 we show entropy colormaps of both models at 
160 ms after bounce. Fluid velocity vectors are superposed, 
providing a snapshot of the flow. By 160 ms after bounce, in 
the nonrotating model s20.nr convection in the high-entropy 
(0(10) ke/baryon) gain layer has developed fully. The shock 
sits at ~ 175 km and is slightly deformed by the onset of the 
SASI. Not visible on the scale of this figure is the lepton- 
gradient-driven convective r egion deep inside the P NS, which 
was extensively discussed in lDessart et al.l (l2006a|) . 

The PNS in the rapidly rotating model s20.7r is rota- 
tionally flattened, with unshocked low-entropy inner-core 
pole/equator asymmetry ratios below ^0.5. The shock 
is slightly prolate and has attained an average radius of 
~230 km. The moment-of-inertia-weighted mean period 
of the unshocked (specific entropy ^ < 3 ke) inner core is 
~2.0 ms. Differential rotation between ^20 and 200 km 
is very large, with the angular velocity ft dropping from 
^1600 rad s~^ to a mere ~15 rad s~^ over this radial equa- 
torial interval. Yet, the specific angular momentum j is still 
monotonically and rapidly increasing. It flattens, but does not 
decrease, only at radii greater than ~100 km. This positive 
gradient in j stabilizes the postbounce core against convective 
instability at low latitudes (Fryer & Heger 2000), confining 
overturn to the polar regions and large equatorial radii where 
the j gradient is less steep. 



4. RESULTS: SNAPSHOTS 

In this section, we present our S„ multi-angle transport re- 
sults for steady-state model snapshots at 160 ms after core 
bounce. We diagnose the angle-dependent neutrino radiation 
fields and carry out a comparison between multi-angle and 
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Fig. 3. — Polar plot of the normalized specific intensity 
/iy(r,i?,(^)/max[/iy(r,'j?,(y^)] in model s20.nr at 160 ms after core bounce, at 
selected equatorial radii, and for Ve neutrinos at £i/ = 12.6 MeV. At each ra- 
dius, we normalize the specific intensity by its local maximum. Shown is the 
variation with ^ at fixed 1^ = 7r/2. The graphs are based on a S 15 calcula- 
tion. At r = 30 km, the radiation field is practically isotropic, but is already 
appreciably forward-peaked at the neutrinosphere (r^y = 55 km; optical depth 
r = 2/3) and thereafter smoothly transitions over ~ 200-3 00 km to the free- 
streaming limit. 

MGFLD transport results based on local and global radiation- 
field variables. 



4. 1 . Angular Distributions 

The quintessential problem in treating neutrino radiation 
transport in core-collapse supernova cores is the fact that the 
neutrino transport mean-free path, the average distance a neu- 
trino can travel without experiencing scattering or absorption, 
changes by orders of magnitude from inside to outside. More- 
over, the neutrino transport mean-free-path A^^ varies locally 
strongly with neutrino energy (ex el) and matter density. As 
a consequence, gray transport schemes are problematic, since 
neutrino-energy averages can be defined only locally and the 
mean neutrino energy varies significantly throughout the su- 
pernova core. 

From a more geometric point of view, the radiation field 
in momentum space goes from being completely isotropic 
(net flux ~zero) to being focussed into the radial direc- 
tion ("forward-peaked") in the free- streaming regime. In the 
MGFLD approximation, the mean intensity /j^ is evolved in 
time and the angular information, in particular the informa- 
tion on the degree of forward-peaking, is captured only by 
computing spatial gradients in Jj^ and employing a flux lim- 
iter to interpolate between diffusion and free streaming. 

The S„ Boltzmann solver in VULCAN/2D is able to self- 
consistently handle the transition from isotropic to forward- 
peaked radiation. Figure [3] depicts the angular distribution 
in the azimuthal angle (p (see Fig. [T]) of the normalized spe- 
cific spectral neutrino intensity /^^ for electron neutrinos at 
12.6 MeV. In Fig.O the polar angle ^ is set equal to 7r/2 and 
the (/:)-distribution is given at various radii in the equatorial 
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Fig. 4. — Hammer-type interpolated (smoothed) map projections of the normahzed specific intensity Iu(i^,^)/Ju in model s20.nr at 160 ms after bounce. 
The colormap is logarithmic and in each individual projection is setup to range from max(/iy(i?,(^)//iy) (red) to 10~'^max(/iy(i?,(^)//iy) (black). Shown is the 
specific intensity of Ue, i^e, and 'V^" neutrinos at Su = 12.6 MeV (rows) on the equator (^=90°, measured from the pole) and at radii of 60, 120, and 240 km 
(columns). The Hammer projection is set up in such a way that i? varies in the vertical from 0° (top) to 180° (bottom) and Lp varies horizontally from -180° (left) 
to +180° (right). Grid lines are drawn in i?- and (/^-intervals of 30°. Note (a) the increasing forward-peaking of lu with increasing radius (and decreasing optical 
depth) and (b) that at any given radius Ijy of 'V^" is more forward-peaked than that of the Ue component, which, in turn, is always more forward-peaked than 
the Ue component. This fact is a consequence of a transport mean-free path that varies with species (and energy; not shown here) and is smallest for the electron 
neutrinos. 

plane of model s20.nr. At 30 km from the center, the radia- 
tion field of z^e's at e^, = 12.6 MeV is nearly isotropic, which 
corresponds to a circle in Fig.O With increasing radius (and, 
of course, decreasing matter density) the transport mean-free 
path at fixed Su increases and the radiation field gradually 
departs from isotropy and becomes more and more forward- 
peaked. We define the neutrinosphere as the surface at which 
the optical depth r^^, given by 

is equal to 2/3. At around this r^,, the neutrinos decouple from 
matter and begin to stream freely. At 160 ms after bounce 
in model s20.nr, the 12.6 MeV Ue neutrinosphere is located 
at r ~55 km. As is obvious from Fig. [3l the radiation field 
at the neutrinosphere is not yet dramatically forward-peaked, 
but becomes so with increasing radius. However, complete 
forward-peaking only obtains at radii > 250-300 km, beyond 
which the angular resolution of our S„ scheme becomes sub- 
optimal, even with n=16. However, calculations with varying 
number of i9 (and, hence, (p) angles reveal that the transition 
from isotropy to moderate and large anisotropy is adequately 
reproduced at small and intermediate radii (out to ~200 km) 
even in the case of Sg. 

For the purpose of displaying and studying the local neu- 
trino radiation field, we provide equal-area Hammer-type map 



projections (lHammedll892|) . Such map projections are new 
to the field of neutrino radiation transport and beautifully re- 
veal the multi-D angular-dependence of the radiation field. 
In Fig. m we present such Hammer projections on the equa- 
tor (spatial = 90°) of model s20.nr at radii of 60, 120, and 
240 km for the three neutrino species included in our sim- 
ulations at Sjy = 12.6 MeV . In each plot, we normalize the 
specific intensity to the mean intensity to set a common scale. 
The colormap is logarithmic and chosen to have regions on 
the sphere with high intensity appear red and regions of low 
intensity appear black. 

For neutrinos on the equator, the momentum- space forward 
direction is (i9 = 90° ^(p = 0). Electron neutrinos generally have 
the shortest transport mean-free path of all species in the core- 
collapse context and decouple from matter at the lowest den- 
sities. The Hammer projection in the top-left corner of Fig.H] 
of the e^y = 12.6 MeV equatorial radiation field at 60 km corre- 
sponds roughly to the blue line graph in Fig.[3l which portrays 
only its variation with (p. At fixed neutrino energy group e^, 
electron anti-neutrinos and ' V^" neutrinos decouple at smaller 
radii. Hence, as Fig. |4l shows, at 60 km, they already manifest 
greater local anisotropy than the z^^s. This trend continues at 
all considered radii in Fig. ID 

In Fig. [51 we again present Hammer projections of the nor- 
malized specific intensity, but this time consider only z^^s. 
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Fig. 5. — Hammer map projections of the interpolated (smoothed) normahzed specific intensity Iu('&,'-p)/Ju at 160 ms postbounce in model s20.nr. The 
projections are set up in identical fashion to Fig. |4] Shown here is the variation of the angular distribution with energy group (columns) and angular position 
(rows) for electron neutrinos. The radius is fixed to 150 km. As expected in the coordinates used for the S^ transport in VULCAN/2D (see Fig.[l), I^, becomes 
forward-peaked into i? = 0° and degenerate in (f along the pole (0 = 0°), forward-peaked into i9 = 45°,Lp = 0° on the diagonal (0 = 45°), and forward-peaked 
'J? = 90° ,(p = 0° on the equator (6 = 90°). The degree of the radiation anisotropy and its variation from forward-peaked at £u = 12.6 MeV to less forward-peaked 
at £u = 37.5 MeV is apparent. 

specific intensity and is defined by eq. (O. In the following, 
we use its normalized variant k^ = Ky/Jj,. 

In spherical symmetry, kj, is diagonal and has a single in- 
dependent component, the Eddington factor ky. For isotropic 
radiation, ky = \/?> and k^^ = diag(l/3, 1/3, 1/3), while in the 
streaming regime, ky = l and k^ = diag( 1,0,0). In the transi- 
tion from isotropy to free streaming, k^ generally varies from 
1/3 to 1, but in special cases, e.g., enhanced radiation perpen- 
dicular to the radial direction, may assume values below 1/3. 
Note that one of the common assumptions of MGFLD is the 
Eddington closure, setting k^ = ^ everywhere. 

In axisymmetry and ignoring velocity-dependent terms, the 
Eddington tensor has four independent components whose in- 
dividual meaning depends on the coordinates chosen^. We 
assume and transform to spherical coordinates for our discus- 
sion, since they make the interpretation of the components 
most straightforward. 

In Fig. [6l we present radial profiles of normalized Edding- 
ton tensor components at selected electron-neutrino energies 
Ey in models s20.nr and s20.7r. The nonrotating model can be 
considered nearly spherically symmetric, and, hence, should 
and does exhibit the expected Eddington-f actor systematics. 
At small radii and high densities, where neutrinos and matter 
are in equilibrium, krr = kr^ro = k^p^p = ^ and with increasing 
radius, krr -^ 1 and {k^^r^^k^^} -^ 0. As expected from the 
basic decoupling hierarchy, the value of the Eddington tensor 

^ Off-diagonal com ponents of the Eddington tensor can be related to ra- 
diation shear viscosity dMihalas & Mihalaslll984l) . which we do not consider 
here. 



keep the radius fixed at 150 km, and vary the neutrino en- 
ergy and the angular position on the grid. The bottom row 
of Fig.[5]shows the normalized I^ at the equator (0 = 90°) and 
for the 12.6-MeV and 35.7-MeV u^ energy groups. The center 
and top rows show the same groups at ^ = 45° and at 6> = 0°, 
respectively. From the discussion of Fig. [H we are already 
familiar with the overall radiation field geometry. The trans- 
port mean-free path scales roughly inversely with e^. Hence, 
at any given position in the postbounce supernova core, more 
energetic neutrinos should be locally more isotropically dis- 
tributed in momentum space than less energetic ones. The less 
forward-peaked angular ly distribution of the higher-energy 
neutrinos reflects this. 

The degree of forward-peaking in i9 and (p of the radiation 
field in the quasi- spherically symmetric nonrotating model 
s20.nr is essentially independent of the angular position on 
the grid and the radiation fields at any given radius can be 
transformed into one another by simple rotation. Because of 
the aspherical and oblate distribution of matter in the rotating 
model s20.7r, the forward-peaking is also a function of polar 
angle. Due to its PNS's oblateness (see Fig. [2]), the neutri- 
nos generally decouple at significantly smaller radii near the 
pole than near the equator, in turn leading to more strongly 
forward-peaked radiation fields in the polar than in the equa- 
torial regions (Janka & Moenchm eveiji il989aubi: iWalder et al.i 
120051: lDessart7rai2()06R 12007^ " 

4.2. Eddington Factors 

The radiation-pressure tensor K^, also known as the Ed- 
dington tensor, represents the second angular moment of the 
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components is a strong function of Sj,. Lower-£^ neutrinos de- 
couple at higher densities, and, hence, have Eddington tensor 
components which depart from | at smaller radii than i/^s of 
higher energy. This systematics applies, of course, to i^^s and 
'V^"s as well. The off-diagonal component kr'd is zero in the 
isotropic region, does not exhibit clear systematics, and stays 
an order-of-magnitude smaller than the diagonal components 
for all Su and species. 

The rotating model s20.7r has a postshock configuration that 
is far from spherically symmetric (Fig. [2]). We present in 
Fig. [6] separate plots for its Eddington tensor components in 
regions near the pole and near the equator. In the polar re- 
gions and at small radii (r<100 km), the Eddington tensor 
components show the same qualitative behavior as in model 
s20.nr. At larger radii, however, the systematics are reversed 
and lower-£j, electron neutrinos have more isotropic radiation 
fields (smaller krr) than their higher-Sj, counterparts. An- 
alyzing their radiation fields and matter coupling in detail, 
we find that this surprising feature is a consequence of elec- 
tron capture and the polar compactness (large density gradi- 
ent due to rotation) of the supernova core. Electron capture 
near the shock leads to isotropic neutrino emission that can 
locally isotropize the radiation field in semi-transparent re- 
gions. With decreasing density and temperature, the mean en- 
ergy of neutrinos emitted by capture processes shifts to lower 
Siy. This leads to greater local isotropization of lower-£j, neu- 
trinos, which in turn is reflected in the more slowly increas- 
ing krr of these neutrinos. This interpretation is confirmed by 
the fact that we do not find any such feature in the Edding- 
ton tensor components of the ' V^" neutrinos that are not pro- 
duced in capture processes. We also do not observe significant 
isotropization in the Ue radiation fields, since the emission of 
UeS by positron capture on neutrons is weaker due to the lower 
positron abundance. 

In regions of model s20.7r near the equator where the PNS 
is most extended, the neutrino radiation fields stay isotropic 
to large radii and decouple from matter only slowly with ra- 
dius. Since the matter densities in the equatorial plane stay 
roughly a factor of four larger than in the polar regions, the 
cross-over feature in {krr^kr^r^ , k^^ } does not appear and these 
components follow the standard decoupling hierarchy. Inter- 
estingly, and different from in the nonrotating model, k^^ and 



k^^ show quantitatively distinct variation with e^, the latter 
exhibiting significantly less variation with Su at any given ra- 
dius. The interpretation of this observation is not straight- 
forward, but we suggest that it can be attributed to the fact 
that in model s20.7r the radiation field at any given point on 
the equator of the rotationally-flattened core and for any e^, 
and neutrino species varies locally less in the 13 direction than 
in the (p direction. This, in combination with the fact that 
on the equator the radiation field asymptotically peaks into 
the ("3 = 0,(p = 0) direction, results on average in smaller k^^ 
with less spread in energy than exhibited by k^^. The off- 
diagonal component k^^ (not shown in Fig. [6]) vanishes in 
krr = k^'d = k^ip = I regions, but can become relatively large 
at greater radii (up to ~0.2 in magnitude; increasing with e^, 
and radius) and flips sign at the equator. The interpretation of 
kr'd is not straightforward, since its magnitude depends on the 
choice of coordinates. We do not attempt to study it, nor its 
implications for neutrino shear viscosity, in any detail. 

4.3. Global Radiation Field Diagnostics: Luminosities, 
Spectra, Flux Factors, and Neutrino Energy Deposition 

So far we have studied aspects of neutrino transport inac- 
cessible to MGFLD. We now go on to discuss radiation field 
diagnostics that facilitate a S„-MGFLD comparison. F or fur- 
ther reference and compariso n with previous studies (iJankal 
il992; Messer et al. 1998; Burrows et al.ll200Ql) . we define the 
neutrino luminosity per species Lj^. at spherical radius r. 
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where iv is the spectral radial neutrino flux in species vt at en- 
ergy £jy. duo is the spatial solid-angle element, duo = 27r sinOdO 
in axisymmetry. Furthermore, we define the mean inverse flux 
factor (1/F^.), 

1 \ cJdSj,E(s^,Ui) 
V^.j IdSyFr^ey.Vi) ' 

where Eie^^v^ - \i\c~^J{ey^v^ is the spectral neutrino en- 
ergy density, and the neutrino RMS energies are 
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Fig. 7. — S«-MGFLD comparison for the nonrotating model s20.nr at 
160 ms after bounce. All S^ results where obtained with a 16 i?-angle cal- 
culation. See text for details and discussion. Top panel: Neutrino luminos- 
ity as a function of radius and broken down into the three neutrino species 
considered. The 'V^" neutrinos dominate in luminosity and their luminos- 
ity profiles are scaled by a factor of 1/4 to preserve the overall scale of the 
plot. Center panel: Angle-averaged energy-mean inverse neutrino flux factor 
profiles. Bottom panel: RMS neutrino energy profiles. 

The above three quantities are particularly useful diagnos- 
tics, since the e^, -averaged energy deposition rate by charged- 
current absorption of Ve an d Ve on neutrons and protons scales 
linearly with their product (iMesser et al.lll998h . 



4.3.1. Model slO.nr 

In Fig.|71 we plot neutrino luminosities Lj,., angle- averaged 
mean inverse flux factors, and the angle- averaged £^rms 
for the postbounce snapshot at 160 ms of the nonrotating 
model s20.nr. The asymptotic total luminosity at this time 
is ~150 B s~^ and is already dominated by the thermally- 
produced 'V^"s that cool the PNS, but contribute little to 
the heating in the gain region, since they cannot take part 
in charged-current absorption processes. In this quasi- 
spherically symmetric model, we define a spherical gain ra- 
dius Tgain as the radial position beyond which net neutrino en- 
ergy deposition occurs. At 160 ms after bounce, rgain — 90 km 
and the gain region extends almost out to the shock at 
~175 km. The MGFLD luminosities in Fig. |7]are system- 
atically lower by ~5% for i/^s, ~3.5% for i/^s, and ~4% for 
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Fig. 8. — Neutrino luminosity spectra extracted at a radius of 500 km 
for i/g, f?e, and 'V^" neutrinos at 160 ms after bounce in model s20.nr. Solid 
lines correspond to Sn results, while dashed lines are obtained using MGFLD. 
The spectra have the canonical shape and the quantitative behavior found 
in n onrotating intermediate-time postbounce supernova calculations (e.g., 
iTho mpson et al. 2003) with the 'V^ "-neutrinos peaking at the highest en- 
ergies, since they decouple from the fluid at the smallest radii. MGFLD and 
S„ spectra agree closely in shape, but MGFLD is overestimating slightly the 
total asymptotic luminosity (cf. Fig|7). 

'V^"s, but qualitatively resemble the S„ luminosity profiles in 
the gain region. At around the shock position, all MGFLD 
luminosities increase by ~5%. This is a due to the combi- 
nation of the artificially spread-out shock (over ~4-5 zones), 
the rapid change of the inverse neutrino mean-free path in the 
spread-out shock and the implementation of the flux limiter in 
VULCAN/2D. Since this MGFLD artefact occurs right at the 
shock, it can have only little influence on the heating in the 
gain region, but leads to somewhat overestimated asymptotic 
luminosities in the MGFLD case. 

The center panel of Fig. [7] shows the e^, -averaged inverse 
flux factors for the three neutrino species in the MGFLD and 
S„ steady-state calculations of model s20.nr. For isotropic 
radiation (1/Fj,.) tends to infinity, while it approaches one 
when the radiation field becomes forward-peaked at low op- 
tical depth. Focussing on the gain region between rgain and 
the shock position, we find that MGFLD yields mean inverse 
flux factors that are up to ~5% larger for i>eS (less for the 
other species) in the inner gain region. At radii > 150 km, 
the MGFLD (l/F^y.) quickly drops to 1 (free streaming), be- 
coming up to 8% lower than the S„ values in the outer gain re- 
gion. We note that ' V^" interact only via neutral-current weak 
interactions, hence, decouple from matter at higher densities 
and temperatures. Next in the decoupling hierarchy are elec- 
tron anti-neutrinos followed by electron neutrinos. Both S„ 
and MGFLD realize this hierarchy at radii below ~150 km, 
beyond which MGFLD rapidly transitions to free streaming 
irrespective of neutrino species. 

The behavior we observe with radius of the luminosity and 
mean-inverse flux f actor agrees with the general findings of 
iMesser et"aD (Il998|) . In particular, we agree with their assess- 
ment that the artificially accelerated transition to free stream- 
ing in MGFLD occurs not at the neutrinospheres (which are 
generally below the gain region), but at relatively large radii 
within which most of the neutrino source is enclosed. 

In the bottom panel of Fig. |7] we present profiles of the RMS 
neutrino energy for all species in MGFLD and S„ snapshots 
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Fig. 9. — Angle-averaged specific neutrino net gain profile in the s20.nr 
model at 160 ms after core bounce. Shown are the MGFLD results, as well 
as results from steady-state S^ calculations with 8, 12, and 16 'j?-angles, cor- 
responding to a total number of angular zones of 40, 92, and 162. The gain 
region extends from ~90 km to the shock position at ~175 km. The three 
different S« resolutions yield net gain profiles that agree very well (relative 
differences below 1% even for Sg). The MGFLD calculation underestimates 
the total net gain in the outer gain region by at most 10% locally and < 5% 
on average. 
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Fig. 10. — 2D colormap of the specific (per gram) net gain distribution in 
model s20.nr at 160 ms after core bounce. The left half of the plot depicts 
the MGFLD result, Sn is shown on the right. The differences between S« 
and MGFLD are marginal at this time in this model and are practically indis- 
cernible by eye. As a consequence of convection in the gain region and the 
onset of the SASI, even this nonrotating model exhibits significant angular 
and radial variations in the neutrino energy deposition not captured by the 
average profiles in Fig. [9] 

of model s20.nr. The corresponding luminosity spectra (ex- 
tracted at 500 1cm) are shown in Fig. [H Both MGFLD and 
S„ capture the energy systematics that is set essentially by the 
matter temperature in the decoupling region. Neutrino species 
that decouple at smaller radii (higher densities and tempera- 
tures) have higher RMS energies and harder spectra than neu- 
trinos decoupling at larger radii. Quantitative differences in 



RMS energies and in the spectra between MGFLD and S„ are 
small, the slightly higher MGFLD spectral luminosities being 
mostly a result of the artificially enhanced MGFLD luminosi- 
ties near and beyond the shock. 

We now conclude our discussion of the 160 ms-postbounce 
snapshot of model s20.nr by considering the instanta- 
neous neutrino energy deposition rates. Figure [9] depicts 
angle-averaged radial profiles of the specific neutrino heat- 
ing/cooling rates in units of erg (g s)~^ . The region of net gain 
extends from ~90 km to the shock radius and the chief contri- 
bution to the heating comes from charged-current i^^ -capture 
processes on protons, exceeding the corresponding i/^-capture 
on neutrons by a factor of two and more in the narrow ra- 
dial interval from 145 to 175 km. MGFLD underestimates 
the specific net gain in the angle-averaged radial profile by 
at most 10% locally and by ~5% on average at radii greater 
than ~110 km. The integral total net gain predicted by Si6 
is 2.13 B s~^ This is only 3% larger than the MGFLD value 
of 2.07 B s~^ We note in passing that Sg overestimates the 
integrated gain rate by at most ~1.6% while S12 agrees with 
S16 to better than --0.3%. 

Figure [TOl depicts the 2D distribution of neutrino heating 
and cooling in the snapshot of model s20.nr considered here. 
Regions of net gain range from green to red, cooling regions 
are blue to black. The colormap demonstrates the somewhat 
misleading character of angle- averaged profiles. While we 
find that there is little spatial angular variation in the neutrino 
radiation field, the neutrino-matter coupling depends strongly 
on angular position, and energy deposition is generally great- 
est in regions of high entropy (cf. FigO. 

4.3.2. Model s20.7T 

As we discussed in the context of the Eddington tensor in 
^4.21 and as may be guessed from the significant rotational 
deformation of the core in model s20.7r (Fig. O, the radia- 
tion field in this model exhibits a strong rotationally-induced 
asymmetry between pole and equator. In Fig.[TTl we present 
2D colormaps of the radial spectral flux component (in erg 
s~^ cm~^ MeV~^) and isoenergy-density contours (47tJu/c in 
erg cm~^ MeV~^) at a representative Sj^ of 12.6 MeV and for 
all species. Numbers for both S„ and MGFLD are compared 
side by side. The global radiation- field anisotropy system- 
atics are qualitatively similar to what was fou nd in the pre- 
vious MGFLD rotating core-collapse study of IWalder et al.l 
(l2005h . At small radii, the radiation field (energy density) fol- 
lows the density distribution and is oblate, but in the snapshot 
at 160 ms after bounce shown in Fig. [TT] has a pole-equator 
ratio of only 1:2. This ratio increases as the PNS cools and 
contracts. The polar compactness of the PNS core leads to 
a decoupling of matter and neutrinos at smaller radii in re- 
gions near the pole, resulting there in greater spectral fluxes 
at higher neutrino energies and in a prolate distribution of neu- 
trino fluxes and isoenergy-density contours. 

The most striking difference between the S^^ and MGFLD 
radiation fields presented in Fig. [TT] is the former's much 
greater prolateness at large radii for all species (and all ener- 
gies, though we show only Sj, = 12.6 MeV). With the MGFLD 
prescription, the prolateness of the flux is muted and does not 
extend to large radii. Though the radiation fields are smoothed 
out at radii > 150 km by MGFLD, the S^^ fluxes and energy 
densities remain prolate through the entire postshock region 
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Fig. 1 1. — Colormaps of the radial spectral flux at Su = 12.6 MeV of Ue (left), Ue (center), and 'V^" (right) neutrinos in the rapidly-rotating model s20.7r 
at 160 ms after bounce. Isoenergy density contours (47tc~^Ju, vertical color legend) are superposed. The left half of each panel displays the MGFLD result - 
Sn is shown in the right half. The radiation fields are oblate in the PNS core and deform to a prolate shape further out. Note that S« predicts a prolateness of 
the radiation field to much greater radii than MGFLD does. The latter leads to nearly spherically symmetric radiation fields at radii greater than > 150-200 km 
independent of neutrino species. 



and beyond. At radii outsi de ~200 km, the typical striping 
pattern of S„ (ICastorll2QQ4|) becomes visible, though not yet 
dominant. 

In Fig. [121 we plot line profiles of the polar and equato- 
rial "luminosities" (Airr^Fr) of each neutrino species. Pro- 
files obtained with S„ and MGFLD are shown. The asymp- 
totic "luminosities" obtained with S„ have pole-to-equator ra- 
tios of 2.2 {ue). 1.8 {ve), and 2.4 ("i^^"). MGFLD smoothes 
out these large asymmetries, yielding higher equatorial and 
significantly lower polar "luminosities" at radii greater than 
~ 100 km. This is consistent with the more qualitative findings 
based on Fig. [Ill We note that the MGFLD variant of VUL- 
CAN/2D still conserves total flux and energy. For the S„ cal- 
culation, we find total asymptotic luminosities of 21.1 B s~^ 
for Ve neutrinos (MGFLD: 20.4 B s"^), 22.7 B s"^ for Ue neu- 
trinos (MGFLD: 22.6 B s"^), and 53.0 B s"^ for "i^^" neutrinos 
(MGFLD: 52.3 B s"i). Hence, S„ and MGFLD total luminosi- 
ties per species agree very well (and differ at most by ~3.5% 
in the Ve case), while their flux distributions disagree signifi- 
cantly. 

Figure [131 depicting polar and equatorial "luminosity" 
spectra (Aur^FriSv)), reveals that in the S„ calculation (polar: 
black graphs, equatorial: red graphs) the neutrino radiation 
emerging from the PNS and postshock environments through 
the polar region not only has greater fluence, but also a signif- 
icantly different and - in the Ve case - a significantly harder 
spectrum. Ue neutrinos decouple at the largest radii. Their 
"luminosity" spectrum observed by a polar observer peaks at 
£jj ~9.5 MeV, while for an observer near the equator it peaks 
at ~7.5 MeV. Both Ve and "i^^" neutrinos (which decouple 
further in) exhibit a smaller variation in peak energy from pole 
to equator. The MGFLD calculation, on the other hand, shows 
much smaller variations in neutrino energy and flux between 
pole and equator (green and blue graphs, respectively). We 
note in passing that the emerging neutrino spectra of model 
s20.7r are systematically softer by up to ~10% in each species 
than those of the nonrotating model s20.nr presented in Fig.[8l 
This is a direct consequence of the rotationally-induced lower 
overall compactness of the PNS in model s20.7r. 

The RMS neutrino energies in model s20.7r show the same 
overall qualitative behavior and decoupling hierarchy dis- 
cussed in the context of model s20.nr. Hence, we do not 



show them here, but rather state quantitative results. They 
do, of course, trace the strong pole-equator asymmetry that 
we observe in the radiation field. The RMS energies in the 
160 ms S„ snapshot are 1 1.5 MeV (pole) and 10.5 MeV (equa- 
tor) for Ve, 16.7 MeV (pole) and 15.2 MeV (equator) for Ve, 
and 25.8 MeV (pole) and 24.8 MeV (equator) for "z/^" neu- 
trinos. The MGFLD values converge at pole and equator to 
10.5 MeV {ve\ 15.2 MeV (i/,), and 25.0 MeV ("i/^"). 

In Fig. [141 we plot polar and equatorial mean inverse flux 
factor profiles for Ve and Ve neutrinos in our steady-state snap- 
shot for model s20.7r. Results from MGFLD and Si6 runs are 
shown. A free- streaming radiation field has an inverse flux 
factor of one. Due to the steeper density gradient in polar re- 
gions, neutrinos decouple from matter at smaller radii than at 
the equator. While MGFLD must handle the decoupling and 
increased forward-peaking of the radiation field via the flux 
limiter, S„ can track it self-consistently. For Ve neutrinos and 
along the poles, S„ predicts significantly greater mean inverse 
flux factors with shallower slopes than MGFLD, indicating a 
more gradual transition to free streaming than predicted by 
the flux limiter. In the radial interval of ~60-100 km, the rel- 
ative difference is ~ 12-19%, decreasing to ~6-12% out to 
200 km. In equatorial regions, the Vg radiation field is some- 
what more forward-peaked in the S„ calculation at radii below 
~120 km, beyond which MGFLD transitions quickly to free 
streaming while S„ approaches it more gradually, exhibiting 
~6-8% larger mean inverse flux factors in the outer postshock 
region. For Ve neutrinos, the behavior of the mean inverse flux 
factors in polar regions essentially mirrors that observed for 
the Ue^- In equatorial regions, the S„ mean inverse flux fac- 
tor of the Ve^ stays below that using MGFLD out to 165 km, 
beyond which the MGFLD Dg radiation field rapidly transi- 
tions to free streaming. At 180 km, the MGFLD v^ mean 
inverse flux factor is ~1% smaller than that predicted by S„. 
At 220 km, this difference has grown to ~5%. 

Having established the overall neutrino radiation- field char- 
acteristics in the 160-ms postbounce snapshot of model s20.7r, 
we now turn our focus to the neutrino cooling and heating 
rates in this model. We have found little difference in the 
net neutrino heating between S„ and MGFLD variants in the 
160-ms postbounce snapshot of the nonrotating model s20.nr. 
However, based on the differences between S„ and MGFLD 
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Fig. 12. — Radial neutrino "luminosity" profiles (47rr^/v) as seen by 
observers near the pole (solid lines) and near the equator (dashed lines) in 
model s20.7r at 160 ms after bounce. Red graphs correspond to S^ results, 
black graphs depict MGFLD results. Top, center, and bottom panels show L^ 
for Ue, Ve, and 'V/^," respectively. All S« results were obtained with n= 16, 
but for comparison we also plot in the top panel polar profiles that were ob- 
tained with Sg and S12 and find that both S16 and S 12 are very well converged, 
while Sg has troubles at radii greater than ~200 km. However, it agrees very 
well at smaller radii with the higher-resolution S,^ calculations. 

in neutrino fluxes, RMS energies, and flux factors we have 
highlighted in this section, we may expect to find significant 
differences in the neutrino heating rates for model s20.7r. 

Figure[T5]depicts 2D colormaps of the neutrino energy gain 
and loss rate per unit mass (accounting for all energies and 
species), computed for the 160-ms postbounce snapshot of 
model s20.7r using S„ (left panel) and MGFLD (right panel). 
At low latitudes near the equator, S„ and MGFLD agree very 
well to the eye. In regions near the pole, both MGFLD and 
S„ show a pronounced region of net loss at z-coordinates be- 
tween ~40 and ~80 km, beyond which a region of net gain 
(colors light blue and green to red) prevails out to the shock 
position at ~230 km. While the gain region has roughly the 
same physical extent in MGFLD and S„, the latter yields sig- 
nificantly higher energy deposition rates. This is particularly 
the case in the lower gain region at polar angles below ~20° 
and at radii between ~80 and 150 km, where the S„ gain rate 
is larger by a factor of two and more. The top panel in Fig. [16] 
provides a more quantitative comparison of S„ and MGFLD 
gain/loss rates, since it contrasts average specific gain/loss 
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Fig. 13. — Top: Neutrino "luminosity" spectra (47Tr^Fr(£u)) in S^ (red) 
and MGFLD (black) variants of model s20.7r as seen by observers near the 
pole. Ve spectra have solid lines, Ue spectra are shown in dashed lines, and 
'V^"s have dashed-dotted spectra. The spectra are taken from a S16 calcula- 
tion at a radius of 300 km at 160 ms after core bounce. Bottom: "Luminos- 
ity" spectra seen by equatorial observers. 

profiles obtained from polar and equatorial 20° wedges. In the 
polar region, the S„ gain region begins at a radius of ~80 km 
(MGFLD: ~88 km) and the S„ specific gain rate magnitude 
exceeds the MGFLD numbers by a factor of 2.6 at 100 km, 
increasing to 3.2 at 200 km. Near the equator, net energy de- 
position occurs only in a small radial interval of ^90-120 km 
and the MGFLD specific gain rate is larger by 80% at 95 km, 
41% at 100 km, and 26% at 110 km. The net energy loss 
between -- 120-2 10 km (captured by both S„ and MGFLD) 
results from strong electron capture that dominates energy de- 
position by neutrino absorption. 

The observed local differences in neutrino energy deposi- 
tion between S„ and MGFLD are due primarily to the vastly 
different degree to which the two schemes capture the global 
pole-equator asymmetry of the radiation field in the rapidly- 
rotating postbounce supernova core of model s20.7r. S„ 
yields much larger fluxes in the polar direction than MGFLD, 
but predicts lower neutrino fluxes in equatorial regions (cf. 
Fig. [T2I). Differences in the radial mean-inverse flux factors 
and RMS energies are much smaller, and, hence, are of only 
secondary importance. The S„ steady-state snapshot yields 
an integrated gain rate of 1.603 B s~^ while MGFLD predicts 
1.637 B s~^ for the s20.7r snapshot under consideration. This 
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Fig. 14. — Mean inverse flux factors in model s20.7r at 160 ms after bounce 
in polar regions (solid lines) and equatorial regions (dashed lines). Shown 
are profiles for Ve neutrinos obtained with Sn (red) and MGFLD (black), as 
well as profiles for Ue neutrinos (S« blue, MGFLD green). Sn and MGFLD 
graphs agree well inside ~50 km at the pole and inside ~80 km in equatorial 
regions. For Ve neutrinos, S^ yields systematically larger mean inverse flux 
factors in polar and equatorial regions. For Ue, however, S^ predicts larger 
mean inverse flux factors in polar regions, yet transitions slightly faster than 
MGFLD to free streaming in equatorial regions. 

corresponds to ~2.1% more energy deposition per unit time 
in the MGFLD calculation. Given the above discussion, the 
reader may be surprised by these numbers. The explanation 
consists of two factors. Owing to rotation, the amount of mass 
per unit volume (i.e., the rest-mass density) is higher at any 
given equatorial radius than at the same radius in the polar di- 
rection. Plotting the neutrino gain/loss rate per unit volume 
instead of per unit gram, the bottom panel of Fig. [16] clearly 
shows the rotation-induced enhancement of the energy depo- 
sition (per unit volume) near the equator and the larger gain 
rate per unit volume predicted by MGFLD at small to inter- 
mediate radii. The second factor is the simple fact that the 
volume of the equatorial gain regions is much larger than that 
of the polar gain regions. 

As we shall discuss in the following section, the large local 
differences in neutrino heating between the S„ and MGFLD 
snapshots have a dynamical consequence for the rapidly ro- 
tating model and lead to a significant polar expansion of the 
shock in the S„ postbounce evolution calculation. 

5. RESULTS: EVOLUTION CALCULATIONS 

In order to study differences between S„ and MGFLD in 
a time-dependent postbounce setting, we follow our relaxed 
160-ms S„ models in fully coupled radiation-hydrodynamics 
fashion for ~340 ms (model s20.nr) and 390 ms (model 
s20.7r) of postbounce time. In parallel with the S„ runs, we 
continue their MGFLD counterparts for the same time span. 

5.1. Model s20.nr 

Since we begin the MGFLD and S„ calculations from an 
identical hydrodynamic configuration at 160 ms after bounce, 
any qualitative or quantitative differences in their evolutions 
must ultimately be due to differences in the neutrino heating 
and cooling between S„ and MGFLD. 

In the left panel of Fig. [T71 we display the time evolu- 
tion of the integral neutrino energy deposition (net gain) in 
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Fig. 15. — Colormaps of energy- and species -integrated specific neutrino 
energy deposition and loss rates in the rotating model s20.7r at 160 ms after 
core bounce (in units of erg s~^ g~^). The left section of the plot depicts the 
MGFLD result and the right shows the result of the Sn calculation. Note the 
distinctively enlarged polar gain regions and greater specific gain of the S^ 
result compared to the MGFLD calculation. This is in part a consequence of 
the larger polar neutrino fluxes and overall greater flux asymmetry in the Sr 
model (see Fig.FTll. A feature prevalent in both S^ and MGFLD versions of 
this rapidly rotating core is an extended loss region between the shock and the 
small gain region at low latitudes (cf. Fig. [16). The material in the loss region 
is still proton rich (Fe>0.4) and efficiently captures electrons as it advects in, 
radiating away a significant flux of neutrinos (see, e.g., the increase in the 
equatorial "luminosity" between 120 and 150 km in the Sn variant of this 
model, visible in the top panel of Fig.FT2l. Note that both MGFLD and S^ 
exhibit a very small artefact (lower gain/loss) at the symmetry axis associated 
with imperfect numerics/regularization. 

the gain region of model s20.nr. The net gain systematically 
declines at early postbounce times, due (a) to the declining 
neutrino luminosity and (b) to the rapid settling of accreting 
mat erial into the net loss r egion near the PNS core (cf. Fig. 
7 of iMarek & Jankall2007h . At later times, SASI-modulated 
convection increases the dwell time of accreting outer core 
material in the gain layer and the slope of the net gain evo- 
lution flattens. Both S„ and MGFLD track these systematics 
without qualitative difference. The S„ calculation predicts on 
average ~5-10% higher net gain in the postbounce interval 
from ~160 ms to ~220 ms. Between ~220 ms and ~280 ms, 
MGFLD and S„ net gain rates agree to within a few percent. 
Towards the end of this interval, the net gain of the S„ calcu- 
lation grows and settles at values that are on average 20-30% 
higher than those of the MGFLD run. This trend is confirmed 
by the right panel of Fig. [T71 which portrays the heating effi- 
ciency, defined as the ratio of net gain to the sum of Ve and Ve 
luminosities. 

The left panel of Fig. [18] depicts the temporal evolution of 
the Ve "luminosities" (Anr^Fr) as seen by observers situated 
at 250 km along the north pole and south pole as well as in 
the equatorial plane of models s20.nr and s20.7r. Here we 
focus on model s20.nr and note for the S„ variant that north 
pole (thin solid black lines) and south pole (thin solid green 
lines) "luminosities" agree (on average) in magnitude, but ex- 
hibit oscillations about their temporal average that are roughly 
out of phase by half a cycle. The MGFLD calculation (thin 
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Fig. 16. — Top: Averaged specific radial neutrino gain and loss profiles 
in model s20.7r at 160 ms after core bounce. Shown are results from the 
S„ (red) and MGFLD (black) calculations. Both polar and equatorial radial 
profiles are obtained by averaging over 20° wedges. As is already clear from 
Fig. [15] Sn yields significantly greater polar specific neutrino energy gain than 
MGFLD. The S^ gain region extends further in by ~10 km and the gain is 
more than a factor of two larger in the interval fro m ^ 90 to 200 km. Given the 
larger flux asymmetry in the S^ calculation (Fig.[TT), less neutrino flux is go- 
ing through regions of low latitude, resulting in the lower specific gain at low 
latitudes predicted by S^. Bottom: Neutrino gain density (density-weighted 
specific gain). Due to rapid rotation higher densities obtain out to larger radii 
at low latitudes. This results in a partial reversal of the picture presented by 
the upper panel; weighted by density, the neutrino gain (now per unit vol- 
ume) in the equatorial wedge becomes comparable to that near the poles. 
Furthermore, equatorial regions, since they subtend the largest solid angles, 
contribute most to the volume integral. The integral numbers for the net gain 
in the polar wedge (counting both poles) for S^ (MGFLD) are 0.17 B s~^ 
(0.047 B s~^) and in the equatorial wedge are 0.35 B s~^ (0.47 B s~^). The 
total integrated net gain is 1.603 B s~^ and 1.637 B s"^ for S„ and MGFLD, 
respectively. These numbers are surprisingly close given the large qualitative 
and quantitative local differences in the neutrino gain distribution. 



dashed lines), on the other hand, does exhibit some short- 
period "luminosity" variations, yet shows no appreciable dif- 
ference between poles and equator. 

The time at which S„ begins to yield systematically larger 
neutrino heating rates (Fig. [17]) coincides with the growth 
of the SASI-related shock excursions to large amplitudes 
(Fig. [19]). This suggests that the increased heating is related 
at least in part to the S„ variant's ability to better capt ure ra - 
diation field asymmetries (see also the discussion in ^4.3.21) , 
induced at late times by the rapidly varying shock and post- 



shock hydrodynamics in this model. Other factors that con- 
tribute to the increased heating in the S„ calculation are the 
higher RMS neutrino energies (by ~5%; shown in the right 
panel of Fig. [18]) and the more gradual transition of the S„ 
neutrino radiation field to free streaming in the postshock re- 
gion (see ^4.3.1t . 

Figure [20l contrasts S„ and MGFLD simulations of model 
s20.nr by means of colormaps depicting the specific entropy 
distributions in the two variants. To visualize the hydrody- 
namic flow, we superpose fluid velocity vectors. Each panel 
of this figure corresponds to a specific postbounce time and 
each panel's left-hand- side depicts the state of the MGFLD 
calculation, while the right-hand- side depicts the correspond- 
ing S„ calculation. The figure covers a postbounce interval 
from 160 ms (top left) to 500 ms (bottom right). At the be- 
ginning of the runs, the SASLdriven deviation from spheric- 
ity of the stalled shock is mild, but grows with time, show- 
ing i = I excursio ns now generally recognized as character- 
istic of the SASI^Q dScheck et al.ll20()8l: iMarek & Jankal l2()()7l: 
iBruenn et aDl2006l: iBurrows et al.ll2007ch" 

As expected from the discussion of the s20.nr 160-ms post- 
bounce steady-state snapshot in ^ S„ and MGFLD variants 
of this model do not differ significantly in the early SASI 
phase. However, at later SASI stages, in particular at post- 
bounce times > 300-350 ms, the simulations diverge, show- 
ing different local qualitative and quantitative behavior within 
the overall SASI theme. This is also reflected in Fig. [19] 
which depicts the evolution of the average shock radius, as 
well as the shock radii along north pole and south pole. The 
shock positions in the S„ and MGFLD simulations remain 
close and the SASI stays practically in phase (right panel of 
Fig. [19]) until ~ 350 ms after bounce. Only then do they be- 
gin to show significant departures from each other. The SASI 
in the S„ calculation appears more pronounced at later times, 
exhibiting larger local (in time) shock excursions. Yet, quite 
surprisingly, given the significant increase in neutrino energy 
deposition, the S„ calculation does not exhibit any increase in 
the average shock radius, nor does it appear to be any closer 
to explosion than its MGFLD counterpart. 



5.2. Model s20.7r 

The diagnosis of the radiation-hydrodynamic evolution of 
the rapidly spinning model s20.7r is less straightforward than 
for the nonrotating model s20.nr. As discussed in ^4.3.21 ro- 
tation creates a global pole-equator asymmetry in the hydro- 
dynamics of this model. MGFLD and S„ track the effect of 
globally asymmetric matter distributions on the neutrino ra- 
diation field to different degrees. In the steady-state snapshot 
at 160 ms, S„ predicts stronger neutrino heating in polar re- 
gions, yet weaker heating in the higher-density, larger- volume 
equatorial regions. 

The polar, equatorial, and angle-averaged shock positions 
portrayed by Fig. [19] show that the hydrodynamics responds 



^^ At least in detailed 2D models. 'Iwakami et aU I2OO80 carried out an 
exploratory 3D numerical study with nonrotating progenitors that suggests 
that in the 3D case the £=1 dominance still obtains, yet reaches smaller rela- 
tive amplitudes, since not only higher i modes, but also m modes, may now 
contain power. However, Yamasaki & Foglizzo (2008), who performed a per- 
turbative study without symmetry constraints, argued that in the 3D case with 
rotation, a dominant m = 1 (m = 2) mode is likely to emerge in the case of slow 
(rapid) rotation. 
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Fig. 17. — Left: Evolution of the total neutrino net gain rate as a function of postbounce time in the S^ and MGFLD variants of models s20.nr and s20.7r. At 
postbounce times before ~300 ms in model s20.nr, S^ yields a net gain rate that is larger by (on average) ~ 10-15% than that predicted by MGFLD. As the SASI 
becomes more pronounced at postbounce times > 300 ms, the S« net gain begins to more significantly exceed that of MGFLD, averaging out at ~20-30% larger 
values than the MGFLD net gain rate. In model s20.7r, S« and MGFLD net gain rates stay very close in the first ~30 ms of evolution, yet depart when the S^ 
variant approaches its new dynamical equilibrium (see Fig. [19} and provides for a larger gain region (mass and volume). This leads to a net gain rate that is larger 
by ~20-25% (on average in the postbounce interval from 200-350 ms). At later times, the MGFLD calculation, approaching the Sn variant's postshock extent 
(Fig. [19) produces larger net gain rates due to its larger equatorial neutrino fluxes at similar hydrodynamic configuration. Right: Heating efficiency evolution 
in the two models with their S^ and MGFLD variants. We define the heating efficiency as the ratio of total neutrino net gain rate and the sum of electron and 
anti-electron neutrino luminosities. 
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Fig. 18. — Left: Ve "luminosities" (47Tr^Fr) as a function of postbounce time as seen by observers located at a spherical radius of 250 km along the north 
pole (black lines), south pole (green lines), and in the equatorial plane (red lines) in S^ (solid lines) and MGFLD (dashed lines) variants of model s20.nr (thin 
lines) and s20.7r (thick lines). Note that the south pole, north pole, and equator MGFLD "luminosities" in model s20.nr (thin dashed lines) are very similar. Their 
lines are indistinguishable. The same holds for the south and north pole MGFLD "luminosities" in model s20.7r (thick black and green dashed lines). Right: 
Angle-averaged RMS energies of i^e (solid lines) and Ue (dashed lines) neutrinos as a function of postbounce time in the S^ and MGFLD simulations the two 
models. S« predicts systematically higher RMS neutrino energies in both models. 



immediately to the increased polar heating in the S„ calcula- 
tion by a pronounced expansion of the shock along the poles. 
This expansion lasts for ~40 ms, after which the shock has ex- 
panded by ~20% from ~230 km to ~275 km on both poles. It 
stagnates at this radius and subsequently contracts again when 
feedback of the hydrodynamics to the neutrino microphysics 
leads to increased cooling (cf. the increased polar neutrino 
emission shown in Fig.fTSl). The increased postshock volume 
also results in a larger gain region and increased (compared 
to MGFLD) total neutrino energy deposition and heating effi- 
ciency. However, this increased heating is not able to sustain 
the large postshock volume. The shock slowly recontracts in 
the postbounce interval from '^250 ms to ~380 ms and even- 
tually settles at radii similar to those obtained by the MGFLD 



shock. 

In Fig. 21, we present a sequence of 2D entropy colormaps 
with superposed velocity vectors, portraying the postbounce 
evolution of model s20.7r from 160 ms on. The rapid rota- 
tion in this model not only partially stabilizes convection, but 
also weakens and delays the growth of the characteristic £=l 
SASI^^ Since larger shock radii are associated_with an_Jn- 
creased gro^yth rat e of the SASI (e.g., Fogliz zoet al.ll2QQ7l: 
IScheck et al.ll2008h . the S„ variant begins to develop periodic 
shock excursions along the symmetry axis at much earlier 
times than the MGFLD simulation (Fig. [T9l) . However, at 



^^ But see the work of lYamasaki & Foglizzol (120081) . who find via pertur- 
bative analysis that in 3D, rotation enhances the development of azimuthal 
m = 1 and m = 2 SASI-related spiral structure. 
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Fig. 19. — Left: Average shock radii as a function of postbounce time in S^ (red) and MGFLD (black) variants of the nonrotating model s20.nr. Also shown 
are the overall average shock radius, the average of south-pole and north-pole shock radii, and the equatorial shock radius for the rapidly spinning model s20.7r, 
again for S« (blue) and MGFLD (green). In model s20.nr, MGFLD and S^ show little quantitative deviation from each other. In the s20.7r evolution, however, 
a significant increase in the various shock radii is noticable right at the beginning of the time-dependent S^ calculation. At later times MGFLD catches up and 
the average shock radii approach each other. The S^ variant exhibits larger variations, indicating stronger SASI-like shock excursions. Right: Evolution of the 
north-pole (positive) and south-pole (negative) shock radii for the S^ and MGFLD variants of the two models. Since the lowest-order and dominant mode of 
the 2D SASI is the i-\ polar sloshing mode, the polar shock radii are good indicators of its strength and periodicity. Note the initial suppression, but late-time 
development of SASI-like polar shock excursions in the rotating model. 

times later than ~400 ms, the MGFLD model picks up the 
large-amplitude SASI as well and both calculations exhibit 
large-scale radial shock excursions beyond ~400 km along 
the pole (Fig. [191). The average shock radius increases in both 
calculations in this late postbounce phase. We observe nei- 
ther such large shock excursions nor a systematic late-time 
increase of the average shock radius in the nonrotating model. 
The observed behavior is most likely due to the rapid rotation 
and the resulting rarefaction of the polar regions that reduces, 
in particular at late times, the ram pressure of accretion and 
allows for the more pronounced SASI. 



In the left panel of Fig. [Tsl we contrast the Ve "luminosi- 
ties" (47rr^Fr) seen by observers located at a radius of 250 km 
above the north pole, the south pole, and in the equatorial 
plane of model s20.7r. The MGFLD variant predicts a pole- 
equator flux asymmetry of < 10% that is roughly constant with 
time. The S„ calculation yields a very different picture. Polar 
and equatorial "luminosities" at 250 km (i.e., near the shock) 
are vastly different (cf. ^4.3.21) . Over time, the equatorial "lu- 
minosity" decreases while the "luminosity" along the poles 
is enhanced. At ~200 ms, polar and equatorial "luminosities" 
differ by a factor of ~3. By ~500 ms, this factor has grown to 
4. In addition, the S„ simulation shows SASI-induced varia- 
tions in north and south-pole "luminosities" that grow to ~3- 
5% at late times and are not tracked in the MGFLD variant. 
These variations are akin those reported for the nonrotating 
model s20.nr, yet have longer periods, since the large shock 
excursions in model s20.7r occur on longer timescales. 

As in the nonrotating model, we also find in model s20.7r 
that S„ yields systematically higher RMS neutrino energies 
for all species and at all times. However, as shown in the 
right panel of Fig. [181 the angle- averaged RMS energies do 
not exhibit a significant increase in the time interval covered 
by our simulations. This, again, is due to rapid rotation which 
slows down the PNS's contraction. Not shown in Fig.[T8l but 
present in the S„ variant throughout its postbounce evolution, 
are ~ 10-20% (roughly constant in time and independent of 
species) higher RMS energies for neutrinos emitted from po- 



lar regions compared to those emitted from the PNS equator. 
This is consistent with our analysis of the neutrino spectra 
and RMS neutrino energies for the 160 ms postbounce steady- 
state snapshot presented in ^4.3.21 

We end our postbounce simulations of model s20.7r with S„ 
and MGFLD at 550 ms after bounce. Though within roughly 
the same qualitative picture, the two approaches to neu- 
trino transport yield appreciable differences in the postbounce 
radiation-hydrodynamics evolutions. Importantly, and in con- 
trast to our findings for the nonrotating model, S^^ in model 
s20.7r does not lead to systematically higher integral neu- 
trino energy deposition, and at late postbounce times, shows 
a volume-integrated heating rate that is even ~30% lower (on 
average) than in its MGFLD counterpart. 

6. SUMMARY AND DISCUSSION 



Using the code VULCAN/2D (JLivneeLaLl [20041: 
Burro ws et al.l l2007d: iLivne et aP l2007h . we perform long- 
term full-2D multi-angle, multi-group neutrino radiation- 
hydrodynamic calculations in the core-collapse supernova 
context. Based on postbounce hydrodynamic configurations 
from MGFLD simulations, we first compute 2D angle- 
dependent (S„) steady-state solutions for models without 
precollapse rotation and with rapid rotation (fto = ir rad s~^). 
From these snapshots, we numerically follow the radiation- 
hydrodynamics evolution with S„ neutrino transport, tracking 
the nonrotating model to 500 ms and the rotating model to 
550 ms after bounce. 

Done for the first time in 2D, we investigate in detail the 
angle-dependent specific intensities and neutrino radiation 
fields. We compute angular moments of the specific inten- 
sity, including the Eddington tensor, and introduce Hammer- 
type map projections to visualize the angle dependence of the 
specific intensity. These we employ to demonstrate the decou- 
pling sy stematics of the neutrinos and the gradual transition to 
free- streaming of the radiation fields with decreasing optical 
depth. 
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Fig. 20. — 2D entropy colormaps portraying the postbounce evolution of model s20.nr between 160 ms (top-left panel) and 500 ms (bottom-right panel) after 
core bounce. Fluid-velocity vectors are superposed to provide an impression of the flow. Each panel's left-hand side corresponds to the MGFLD calculation 
and each panel's right-hand side shows the S^ result. The time of each panel is given relative to the time of core bounce. The sequence of panels portrays the 
canonical development of the SASI in the nonrotating axisymmetric context. S^ and MGFLD evolution agree very well in the early SASI phases, but deviate in 
detail at later times, while still exhibiting the same overall SASI dynamics. 

We compare our S„ simulations with MGFLD counterparts. 
We find for both models and at all times that the S„ specific 
intensity distributions transition less rapidly from isotropy 
to free- streaming in the semi-transparent outer postshock re- 
gions. Sn yields mean inverse flux factors and RMS neutrino 
energies in these regions that are ~10% larger than those ob- 
tained with MGFLD. In the context of the neutrino mech- 
anism of core-collapse supernova explosions, differences in 
the net neutrino energy deposition rates between MGFLD and 
multi-angle S^^ transport are of greatest interest. In the quasi- 
spherical early postbounce phase of the nonrotating model, 
we find that S„ predicts a 5-10% greater neutrino energy de- 
position rate than MGFLD. At later times, when the SASI has 
reached large amplitudes and globally deforms the postshock 
region, we find that S„ yields consistently larger (up to 30% on 
average) energy depositions and leads to significantly larger 
temporary shock excursions around average shock radii that 
do not depart much from those in the MGFLD calculation. 



Convection on small and intermediate scales and SASI on 
large scales, are the key agents of the breaking of spherical 
symmetry in nonrotating (or slowly rotating) core-collapse 
supernovae. While we observe no large qualitative differ- 
ences in the growth and dynamical evolutions of convection 
and SASI between nonrotating S„ and MGFLD models, we 
find that the imprint of the asymmetric hydrodynamics on 
the neutrino radiation fields is captured with greater detail by 
the multi-angle transport scheme. For the late-time, heavily 
SASI-distorted postbounce core, S„ predicts asymptotic neu- 



trino fluxes that have variations with time and angle of 5-10% 
in magnitude. MGFLD is able to capture the temporal varia- 
tions of the neutrino luminosity, but smoothes out the angular 
flux variations at large radii/low optical depths. 

Rapid rotation leads to large deviations from spherical 
symmetry and a rotationally-deformed PNS emits, by von 
Zeipel's law of gravity darkening, a greater neutrino flux 
along its rotational axis than through its equatorial re- 
gions (Janka & Moenchmever 1989a,b; Kotakeetal. 2 0031: 
IWalder et al.ll2005l: iBuras et al.ir2006al: iDessart et al.1l2006bh . 
We find that both 2D MGFLD and S„ yield similar radia- 
tion fields and pole-equator flux ratios at radii smaller than 
~100 km. At larger radii, the MGFLD radiation fields spheri- 
cize and show little pole-equator asymmetry in their asymp- 
totic variables. S„, on the other hand, captures large pole- 
equator flux ratios of up to 4: 1 at late times and predicts polar 
neutrino spectra that are harder in peak energy (RMS energy) 
than on the equator by up to 30% (10-15%) for Vg neutrinos, 
and somewhat less for the other species. All this results in a 
neutrino energy deposition rate per unit mass in polar regions 
that is locally up to ~2.5-3 times higher when multi-angle 
transport is used. This increased polar neutrino heating has 
a dynamical effect on the postbounce evolution, leading to 
rapid shock expansion in the polar regions and an earlier on- 
set of the (initially) rotationally-weakened SASI. However, at 
late times, the SASI in the MGFLD calculation catches up and 
yields shock excursions of a similar magnitude. 
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Fig. 21. — 2D entropy colormaps portraying the postbounce evolution of the rapidly -rotating model s20.7r between 160 ms (top-left panel) and 550 ms (bottom- 
right panel) after core bounce. Fluid-velocity vectors are superposed to relay an impression of the flow and convey the partial suppression of convective overturn 
in regions of positive specific angular momentum gradient. As in Fig. |20] we plot the MGFLD result on the left-hand side and the S^ result on the right-hand 
side of each panel. Easily discernible is the immediate increase in the polar shock radius in the Sn calculation. This is a direct consequence of the increased polar 
neutrino heating in this variant (Figs. [15] and [16). At intermediate times, S« and MGFLD shock positions grow closer, but later on in the postbounce evolution, 
the Sn variant begins to develop larger top-bottom SASI-like asymmetry and polar shock excursions at earlier time than its MGFLD counterpart. 

of flux limiter and wiU consider different progenitor models. 

The core-coUapse supernova problem is one of many feed- 
backs. Larger heating rates and heating efficiencies than 
found in our models appear to be necessary to break the feed- 
back cycle between neutrino radiation fields and hydrody- 
namics, revive the stalled shock, and unbind the supernova 
envelope - if the neutrino mechanism is to obtain in the way 
presently envisioned. Future work will have to go beyond 
the limitation of axisymmetry and must address in detail the 
entire ensemble of possible factors relevant in the supernova 
problem, including, but not limited to, 3D dynamics, multi- 
angle neutrino transport with velocity dependence and inelas- 
tic Ve-e' scattering, progenitor structure, rotational configu- 
ration, magnetohydrodynamics, convection, the SASI, PNS 
g-modes, general relativity, the nuclear EOS, and neutrino- 
matter interactions. 



In summary, our results show that 2D multi-angle neutrino 
transport manifests interesting differences with 2D MGFLD 
when addressing local and global radiation field asymmetries 
associated with rapid rotation and the non-linear SASI at late 
postbounce times. In addition, multi-angle transport results 
in enhanced neutrino energy deposition. The latter is most 
significant in the polar regions of rapidly rotating postbounce 
configurations and affects dynamically the postbounce evolu- 
tion, including the growth of the SASI. However, in the large 
postbounce interval covered by our simulations, the local and 
global differences between multi-angle transport and MGFLD 
calculations do not appear large enough to alter the overall 
simulation outcome. Importantly, the multi-angle models do 
not appear to be closer to explosion than their MGFLD coun- 
terparts. 

Although we neglect velocity-dependent transport terms 
and coupling of neutrino energy bins, we do not expect our 
conclusions to be altered by their inclusion, since they are not 
likely to affect significantly the differences between multi- 
angle transport and MGFLD. Further significant limitations 
of our present study are the neglect of general relativistic and 
MHD effects, the restriction to only one finite-temperature 
nuclear EOS, the limited resolution in momentum- space im- 
posed by the computational cost of multi-angle calculations, 
and the use of two spatial dimensions, plus rotation. In the 
future, we will investigate the dependence of our results (e.g., 
heating rates, radiation- field asymmetries etc.) on the choice 
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